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Castillejo-Dalitz-Dyson (CDD) poles are known to be connected with bound states and resonances. 
We discuss a new type of CDD pole associated with primitives i.e., poles of the P matrix that 
correspond to zeros of the D function on the unitary cut. Low's scattering equation is generalized 
for amplitudes with primitives. The relationship between the CDD poles and the primitives is 
illustrated by a description of the S-wave nucleon-nucleon phase shifts. 
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The poles introduced by Castillejo et al. in Ref. [ij are 
known as the Castillejo-Dalitz-Dyson (CDD) poles. They 
describe ambiguities in solutions to the Low scattering 
equation for amplitudes that satisfy correct analytical 
properties and unitarity. To clarify the physical meaning 
of the CDD poles, Dyson constructed a model \S\ that 
demonstrates the relation of the CDD poles to bound 
states and resonances. 

Some time ago, Jaffe and Low 4] proposed a method 
for identifying exotic multiquark states with primitives 
that appear as poles of the P matrix rather than the 
S matrix. The analysis, performed for scalar mesons [3| 
and nucleon-nucleon scattering [5], revealed primitives, in 
agreement with expectations from the MIT bag model. 
A dynamical model of the P matrix was developed by Si- 
monov [6] and was applied to the description of nucleon- 
nucleon scattering [6l4llj. Other microscopic models of 
nucleon-nucleon forces have also been discussed [T2| - [l^ . 

The recent interest in the problem of nucleon-nucleon 
interactions is connected to new constraints on the equa- 
tion of state (EOS) of nuclear matter, obtained from col- 
lective flow data and subthreshold kaon production in 
heavy-ion collisions [H, [l^ as well as astrophysical ob- 
servations of massive neutron stars [13, • One-boson 
exchange models are in reasonable agreement with the 
laboratory data but predict surprisingly low masses for 
neutron stars in the (3 equilibrium [19,, . Microscopic 
models can provide better insight into the short-range 
dynamics of nucleon-nucleon interactions and the high- 
density EOS. 

In this Brief Report, we clarify the link between the 
CDD poles and the primitives, which can be useful 
for modeling the nucleon-nucleon interactions in the P- 
matrix formalism. 

CDD poles arise when the interaction of particles in- 
cludes intermediate states that are internally different 
from combined-particle states. These are discrete eigen- 
states of the system and basically form other channels in 
the scattering problem. They have also been referred to 
as elementary particle or compound states. 

In the Dyson model, one starts from the scattering of 
two particles, e.g., a nucleon and a pion. A nucleon can 
absorb a pion and can turn into an excited compound 
state Na of mass > y/sQ = m + fi, where m and fi 



are the nucleon and the pion masses. The D function of 
the process can be written as 



i?(s) = A(s) - n(5), 

where, in the relativistic notations. 
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n(s) = -- / '^^{s')-^^ds' 

Jsn S' -S 
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Here, $2(5) = Trk/^/s is the relativistic two-body phase 
space, k is the center-of-mass momentum, is the cou- 
pling constant, and T{s) is the form factor of the NaNir 
vertex. The S matrix has the form 



D{s-iO) 
D{s + iO)' 



(4) 



The poles of A(s) are the CDD poles. They are located 
between the zeros of A(s) (i.e., between and M^_^_i). 

The D function constructed in such a way is the gener- 
alized R function 1] . It has no complex zeros on the first 
Riemann sheet of the complex s plane. It also has no 
zeros on the real half axis {—oo,so), which corresponds 
to bound states, provided D{so) < and sq < M^. 

The simple roots of the equation 



Dis) = 0, 



(5) 



located on the second Riemann sheet below the unitary 
cut, are identified as resonances. In the limit of small g^, 
roots of Eq. ([5]) are localized in the neighborhood of s = 
M^. SRs gives the renormalized resonance mass, while 3s 
determines the decay width Fq, = g^QD{M^)/Ma- 

At the CDD poles 5{s) = mod (tt), the slope of the 
phase is positive. If s-y is a CDD pole, then Eq. ^ gives 
A''^{sj) = and A~^(s^)' < 0. By expanding the D 
function around s = Sj and by using Eq. (|4]), one finds 

6{s-yY = -5D(s^)A-1(st,)' > 0. 

Such behavior is in agreement with the Breit-Wigner for- 
mula according to which isolated resonances drive the 
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phase shift up. In potential scattering, an increasing 
phase is associated with attraction. 

The Dyson model therefore applies to systems with 
attraction where scattering phase shifts increase with in- 
creasing energy. 

The nucleon-nucleon phase shifts, conversely, decrease 
with increasing energy and provide evidence for repul- 
sion. 

In Refs. [iHl] 3D(s) is strictly positive. Softening this 
constraint to 3Z)(s) > allows extension of the Dyson 
model to systems with repulsion: 

Let us consider the scattering of two nucleons through 
compound states, dibaryons, with form factors J^{s) that 
have a simple zero a.t s = Sp > sq — Am? . Consequently, 
'^D{s) ~ (s — Sp)^. Such behavior is presupposed in 
the quark compound bag (QCB) model developed by 
Simonov for the description of nucleon-nucleon interac- 
tions ,6:]. In the QCB model, the separable potential 
generated by the compound 6-quark bags is restricted 
to the bag surfaces. The 5*- wave form factor has the 
form J-{s) — sin(kb)/kb, where k is the center-of-mass 
momentum and b is the effective interaction radius. In 
relativistic notation, the D function of the model has the 
form of Eq. ([ij, while A(s) and the self-energy operator 
n(s) are equivalent to Eqs. © and respectively. 

This analogy allows the techniques developed in 
Ref. [l| to be used to parametrize nucleon-nucleon scat- 



tering amplitudes with functions that have the correct 
analytical properties. 

If ^D{sp) ^ 0, the phase touches, at s = Sp, one of 
the S{s) ~ mod (tt) levels without crossing. However, 
if Eq. ([5]) holds at s = Sp for both the real and the imag- 
inary parts, the phase crosses one of the levels S{s) = 
mod (tt) with a negative slope. This can be verified by ex- 
panding D{s) around s — Sp. By taking Eq. Q and the 
conditions 3D(sp)" > and diD{spy > into account, 
one gets 



s{spY 



2^D{spy ■ 



In potential scattering, a negative slope of the phase shift 
is associated with repulsion. 

The Low scattering equation [2| is modified in the pres- 
ence of primitives. In the systems with n(s) given by 
Eq. ([3]), the scattering amplitude A{s) — e*'''^''^ sin5(s) 
can be represented as follows: 



Ais) 



1'2(g)-F^(s) 

D{s) ■ 



(6) 



This amplitude obeys the generalized Low scattering 
equation 



J 
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$2(S')-^^(S') S' -S 
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which is essentially the dispersion integral representation 
for the inverse denominator function D(s) that accounts 
for the poles, which correspond to the bound states and 
primitives. A{s) and T{s) have simple zeros at s = Sp, so 
the integrand in Eq. ([7]) is a regular function at s' — Sp. 
The bound states and the primitives generate poles at 
Sb < Sq and Sp > sq on the real axis, the coefficient 

2A{spy 

" $2(Sp)J-2(Sp)" 



is positive. 

In the QCB model, the P matrix takes the form 



P = P. 



free 



(8) 



where Pfree is the free P matrix. For the S wave, Pfree = 
kbcot{kb) and P = kbcot{kb + 6{s)). The value of k is 
fixed by the normalization of £'(s). The compound states 
of masses Ma show up as poles of the P matrix. The 
poles of the P matrix split into two groups according to 
their physical nature: 

The first group is related to the bound states and res- 
onances. 



One bound state always exists at D{so) > 0. Addi- 
tional bound states can be generated by compound states 
with masses Ma < ^/so■ 

A characteristic feature of a resonance is the condition 
J^{s) ^ in the neighborhood of s = M^. Equation ([5]) 
can then be used to find a simple pole of the S matrix. 
The roots of Eq. ([S|) that lie on the real half axis (— oo, sq) 
of the second Riemann sheet are virtual states that can 
be related to the compound states also. 

The poles of the second group are related to roots 
of Eq. ^ on the unitary cut in the neighborhood of 
s — Ma- Such poles do not show up as S'-matrix poles 
and cannot be treated as resonances. They are called 
primitives according to Jaffe and Low j4| . If a resonance 
moves from the second Riemann sheet to the unitary cut, 
its singular effect on the S matrix cancels out. As dis- 
tinct from resonances, primitives drive the phase shift 
down and mimic repulsion. 

In the Dyson model, there exist at most one bound 
state and at most one resonance, which are not associated 
with the CDD poles. In the QCB model, there are CDD 
poles related to primitives that do not give rise to bound 
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states or resonances. The neighboring CDD poles squeeze 
masses of compound states that become bound states, 
resonances, or primitives when coupling to the continuum 
is switched on. This is illustrated in Fig. [TJ 
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FIG. 1: (Color online) D-function zeros in the complex s 
plane. The unitary cut is shown by a bold solid line. Com- 
pound states 1, 2, and 3 are eigenstates of the free Hamilto- 
nian. Upon switching on the s channel interaction with the 
continuum, zeros move to new positions shown by crosses. 
Compound states 1, 2, and 3 become bound state, primitive, 
and resonance, respectively. A pair of the CDD poles that 
squeezes compound state 2 of the primitive type is shown by 



The parametrization ensures the existence of the 
deuteron pole for 



7 = -Kn(sd) > 0. 



(12) 



Unphysical zeros of the D function are eliminated by 
constraining Cp. One can easily show that 3I?(s) ~ 3s, 
and that the coefficient of proportionality is positive for 
positive Cp. In this case, D(s) has no zeros for 3s ^ 0. 
The real half axis (— oo, sq) remains to be checked. The 
derivative I?(s)' is positive below the threshold. D(s) 
crosses the real axis at s = s^ < sq. This is the unique 
zero of the D function, provided A(s) has no poles for 
s < So. Let us investigate the zeros of A~^(s). Since 
K-^A-^{sd) = -1/7 < 0, A-i(s)' < and A-i(s) 
has no poles for s < Sp by construction, the condition 
A~^(— oo) < is sufficient to exclude unphysical zeros. 
Finally, Cp satisfies the constraint 



The S-wsLve nucleon-nucleon scattering can be consid- 
ered as an example of dynamics influenced by the CDD 
poles that are connected to primitives. 

The model we discuss is the relativistic extension of 
the QCB model. The P-niatrix formalism is recovered 
with 



J-(s) = - 
■So 



1/4 



sin(fc6) 



kb 



Equation ([3]) for k — 2mh/'iT gives 



Kn(s + iQ) 



kb 



(9) 
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In addition to the interaction through compound 
states, we introduce a contact interaction. This amounts 
to a redefinition of A~^(s) as compared to Eq. ([2]). In 
the case of one CDD pole, the most general expression 
for A~^(s) becomes 



K-^K-\s) = Cp{- 



Sd 



1 
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(11) 



where Cp is a free parameter such that nCpTp — g^n rp — 
8(7r/&)^ is the residue of Pfree in the 5* wave, and Sd = 
MJ is the deuteron pole or the threshold. The D function 
with the contact interaction remains the generalized R 
function. 

In the channel, the phase shift vanishes at Tiab — 
354 MeV. The S matrix according to Eq. (U) is unit in 
two cases: A(s) = oo and ^D{s) = — 3n(s) = 0. The 
poles of A(s) are the CDD poles. At the CDD poles, 
the slope of the phase is positive, which corresponds to 
attraction. The second case 3n(s) — gives repulsion. 

Tiab = 354 MeV is equivalent to fc = 408 MeV. The 
equation jF{s) = gives kb = tt. Thus, we determine 
b = 1.52 fm. Since J'(sp) = if and only if 5Rn(sp) = 0, 
Eq. ([5]) simplifies to A(sp) — 0. A(s) vanishes when s = 
M^. The compound state shows up as the primitive of 
mass Ma = 2y/WT^ = 2047 MeV. 



< Cp < cT"" = 



Sd 



p 



(13) 



In Fig. [2] (a), we show the ^Si phase shift versus the 
proton kinetic energy for Cp = 0.9c™'*'^. This is compared 
to the partial wave analysis data provided by Ref. W^ . 
For thepn system, s — so+2mnTiab, where m„ is neutron 
mass. The CDD pole is located at M = 3203 MeV. The 
pion production threshold is at Tiab = 280 MeV and the 
inelasticity is small up to ~ 350 MeV. 

In Fig.[3](c), we show ^D{s) and 3D(s) as functions of 
Tiab- ^D{s) has one zero below sq, which corresponds to 
the deuteron. The second zero at s > sg with 3£'(s) ^ 
and a negative slope of ^D{s) ensures the crossing of 
the level 6{s) — tt. The third zero corresponds to the 
primitive. 

In the channel, the phase shift vanishes at Tab = 
265 MeV. The same arguments as before give b = 1.76 
fm and M„ = 2006 MeV. 

A^^(s) has the form of Eq. (|lll) . with Sd replaced by 
So — Am?. Near the threshold, kD{s) = —j+l + ikb+. . . 
From the other side, D{s) ^ 1 — iS{k) + ... = 1 — ika + 
where a = 23.56 fm is the scattering length. One has to 
require 



7 = 1 + 



(14) 



D{s) has no zeros for complex values of s. Its derivative 
is positive for real s < so- To avoid unphysical zeros, it 
is sufficient to require 



kD(—oo) 



< kD{so) = 1 - 7 < 0. 



The second inequality is fulfilled, and the first one gives 
Cp < min(c™^'^, c™^'^/(7 — 1)). Since 6 ^ a, this reduces 
to Eq. ([T3|) with Sd replaced by sq. 

In Fig. [5] (b), we show our fit of the pn -^Sq phase shift 
with Cp — 0.9c™'^^ compared to the experimental data 
The CDD pole occurs at M = 2916 MeV. 
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Shown in Fig. [2] (d) are the real and the imaginary 
parts of the D function versus the proton kinetic energy. 

In Figs. [2] (c,d), the real and the imaginary parts of 
the D functions vanish at s = Sp. These are signatures 
of the primitives, along with crossing the levels 5{s) — 
with negative slopes on Figs. [2] (a,b). 

The values of h and Mq, are close to those obtained in 
Refs. 

Benjamins and van Dijk [l3| used the hybrid Lee 



model with one compound state in each channel to de- 
scribe the nucleon-nucleon 5'-wave phase shifts below 
Tiab — 500 MeV and to reproduce parameters related to 
the deuteron and the virtual ^5*0 state. The model does 
not have explicit CDD poles and primitives. However, it 
can be reformulated in terms of the QCB model with one 
CDD pole and two compound states that correspond to 
the primitive and a high- mass resonance 
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FIG. 2: (Color online) ^Si and scattering phase shifts in radians (upper panel) and real and imaginary parts of the D 
functions (lower panel) versus the proton kinetic energy. The solid curves are parametrizations within the relativistic QCB 
model. The experimental phase shifts f2l| are shown by triangles. 



Resonances and primitives do not exist as asymptotic 
states. In Feynman diagrams, propagators of primitives 
l/(s — M^) are multiplied by form factors J-{s). Such 
combinations do not have poles at s = M^. Primitives 
thus do not propagate, though they influence the dynam- 
ics. 

Summarizing, the physical meaning of the CDD poles 
was revisited. In the general case, the neighboring 
CDD poles squeeze masses of compound states related to 
bound states, resonances, or primitives. The primitives 
are P-matrix poles associated with zeros of the D func- 
tion on the unitary cut, which do not show up as poles 
of the S matrix. The Low scattering equation was gen- 
eralized for amplitudes with primitives. The primitive- 
type CDD poles occur in systems with repulsion. In the 



and nucleon-nucleon channels, the CDD poles 
at M = 3203 MeV and M = 2916 MeV are associated 
with the primitives at = 2047 MeV and M„ = 2006 
MeV, respectively. The model we used ensures that the 
D function has the correct analytical properties on the 
first Riemann sheet of the complex s plane and provides 
the partial wave amplitudes that satisfy the generalized 
Low scattering equation. 
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